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Summary 

Singularity conditions that arise during structural optimization 
can seriously degrade the performance of the optimizer. The 
singularities are intrinsic to the formulation of the structural 
optimization problem and are not associated with the method 
of analysis. Certain conditions that give rise to singularities 
have been identified in earlier papers, along with a proposition 
to alleviate the consequences of their presence (refs. 1 to 3). 
These singularities were global in nature, encompassing the 
entire structure. Further examination revealed more complex 
sets of conditions in which singularities occur. Some of these 
singularities are local in nature, being associated with only 
a segment of the structure. Moreover, the likelihood that one 
of these local singularities may arise during an optimization 
procedure can be much greater than that of the global sin- 
gularity identified earlier. This paper provides examples of these 
additional forms of singularities. It also gives a framework 
in which these singularities can be recognized. In particular, 
the singularities can be identified by examination of the stress- 
displacement relations along with the compatibility conditions 
and/or the displacement-stress relations derived in the 
integrated force method of structural analysis. 


Introduction 

Structural optimization methodologies based on mathe- 
matical programming techniques require constraint-gradient 
information (ref. 4). Nonlinear optimization proceeds 
iteratively, and at each iteration a direction vector 0 is 
generated in the design variable space. The formation of this 
direction vector utilizes the gradients. Consider, for example, 
the coefficient matrix [H] used in a method of feasible 
directions, as given by Best, (ref. 4) to find the direction vector 
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where V/ is the gradient of the merit function. 


K = - 0.5 V | [W]V/| (4) 

and k is the number of active constraints used in the computation 
of the direction vector. 

In equation (3), x, represents the ith design variable, and 
gj represents the jth constraint value (e.g., the nondimen- 
sionalized element stress or nodal displacement). Here, we 
consider only sizing design variables (e.g. , the cross-sectional 
areas of truss elements). 

Now, if the rank of the gradient matrix [Vg] is less than 
k, then the derived composite matrix flV^] 7 [Vg] ] will be 
singular. Thus, the coefficient matrix [H] is not defined, and 
any direction generated from equation (1) would be spurious. 
Gradient projection methods that use the matrix \ H] can suffer 
the same consequences, and similar remarks can be made about 
other techniques that use constraint-gradient information in 
this way. Computer codes can react to this situation by pre- 
mature termination (sometimes with indications of multiple 
overflow errors) or divergence. These observations are intrinsic 
to the specification of the constraint formulation in the struc- 
tural optimization problem and are not associated with any 
specific analyzer. 

Clearly, the rank of the gradient matrix [V#] is less than 
or equal to the smaller of the number of active constraints k 
and the number of design variables n. Thus, whenever the 
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number of active constraints is larger than the number of design 
variables, the rank of [Vg\ will be less than k (since n < /c), 
and singularities in [H] will occur. However, as pointed out 
earlier (refs. 1 to 3), whenever the number of active stress 
and displacement constraints is more than the number of dis- 
placement degrees of freedom of the structure (which is often 
significantly less than n ), singularities will also be introduced. 
This condition depends only on the number of active 
constraints, irrespective of where on the structure the active 
constraints appear, and will be referred to as a global 
singularity condition. 

Linear Functional Dependence Among 
Constraints 

Constraint sets formed during structural optimization may 
exhibit functional dependence. Moreover, this functional depen- 
dence may be linear in form. That is, for certain sets of 
constraint functions of the design variable \ (e.g.,{#|(x), 

#2 (x) g K ( x)!)< there exists a set of constants {<*,-[, not all 

zero, such that 

K 

^2 a jXj ( x) = 0 for all x (5) 

./= l 

Note that any individual function gj(x) may be nonlinear 
in the independent variable X- This linear functional dependence, 
in turn, generates linear dependence among constraint gradients 
(which leads to singularities). Differentiating equation (5) with 
respect to x gives 

K 

ctjVgj = 0 (6) 

./= l 

where is the gradient of the function gj with respect to x* 


Constraint Dependence for a Three-Bar Truss 

To illustrate the existence and importance of linear functional 
dependence among constraints, consider a three-bar steel truss 
with a Young’s modulus E of 30 000 ksi and a strength a () of 
20 ksi (fig. 1). For each load condition, the truss can have 
three stress and two displacement constraints. The stresses and 
displacements are functions of the cross-sectional areas of the 
bar elements x and either (or both) type(s) of constraints may 
appear in the active constraint set. It can be shown that, for 
the three-bar truss, the following relationships hold for each 
load condition: 



(o | - (j 2 + <Ji) = 0.0 (7b) 

where Xj(x), the yth nodal displacement, and n,(x), the /th 
element stress, are related to the behavior constraints and are 
functions of the three-component design variable x (the areas 
of the elements). Here, f is the length of the second element 
(see fig. 1). 

Note that these relationships hold for all values of the design 
variable (since x, does not appear explicitly in eqs. (7)). 
Furthermore, linear functional dependence (or independence) 
of a set of behavior variables \o n Xj\ implies the linear 
dependence (or independence) of their associated constraint 
gradients (ref. 5). 

Since we have three equations in five unknown functions 
(<7j, u 2 , a 3’ ^i> an d X 2 ), given any two of the five functions 
(except for the <j 2 , X 2 pair), the other three functions can be 
determined. The observation that any set of more than two 
constraints is linearly functionally dependent corresponds to 



Figure 1 .— Three-bar truss. (FJcments are circled, nodes arc not.) 
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TABLE I -RANKS OF SEVERAL CONSTRAINT-GRADIENT 
MATRICES FOR THE THREE BAR TRUSS 
[Element areas: 1.00, 9.30. and 3.67 in. 2 ] 


Case 

Constraint set 

Rank of 
constraint- 
gradient 
matrix 

Nature of 
singularity 

Number 

Constraints 

1 

5 

1 0 1 , (7-> | ,Xi j 

2 

Global 

2 

3 

|<7| 

2 

Global and local 

3 

2 


I 

Local 







the global singularity condition identified earlier (refs. 1 to 3). 
The linear functional dependence among stresses, also noted 
earlier (refs. 2 and 3), is clearly indicated by equation (7b). 
Use of the knowledge of this global singularity to restrict the 
active constraint set to two constraints is, however, insufficient 
to guarantee linear functional independence. In particular, 
equation (7a. 2) involves only one element stress (<r 2 ) and one 
nodal displacement (X 2 ) that would produce a singularity 
because of their linear functional dependence. We would 
expect, then, that the rank of the constraint-gradient matrix 
formed by (1) taking all five constraints would be two, (2) 
taking the three stress constraints would certainty be two, and 
(3) taking a 2 and X 2 would be one. Table I shows the results 
of performing a singular value decomposition on the three 
constraint-gradient matrices indicated. 

Use of the IFM Equations for Identification of Singularities 

The integrated force method (IFM) is a structural analysis 
tool based on the method of forces through which global 
singularities in structural optimization were identified earlier 
(refs. 1 to 3). It will be used here to identify additional types 
of singularities. 

Let us consider, more generally, an arbitrary truss under 
a single load condition that has n stress degrees of freedom 
and m displacement degrees of freedom. We will designate 
this structure as truss (njn), following the convention in the 
IFM (ref. 6). The following relations can be derived from the 
governing equations used in the IFM (see appendix B): 


a = [B\X 

(8a) 

X=[J\ 0 

(8b) 

0 = [C] 0 

(8c) 


where the matrices [Z?], f/|, and [C] are independent of the 
design variables; have full ranks of m, m, and r, respectively; 
have dimensions of n X m, rn X n, and r x n, respectively; 
and r = n — m. The choice of design variables depends on 
the structural elements chosen for the design, but typically these 
variables are taken as the element cross-sectional areas for truss 
elements, the moments of inertia for flexural elements, and 


the shear area and polar moment of inertia tor elements with 
torsion. 

Equation (8a) is referred to as the set of stress-displacement 
relations and is derived from the force equilibrium equations. 
Equation (8b) describes the displacement-stress relations. 
Equation (8c) is the set of stress compatibility conditions. We 
note here that equations (8b) and (8c) can be used to derive 
equation (8a), and vice-versa. From the stress compatibility 
conditions (eq. (8c)), we can see that the set of all stress con- 
straints is linearly functionally dependent. These equations are 
also consistent with the observation regarding global singu- 
larities made earlier (that the maximum number of stress and 
displacement constraints which can form a set ot linearly 
independent functions is m) (ref. 1). However, more subtle 
cases of singularities can arise by linear functional dependence 
among subsets of stress constraints (see eq. (8c)) or stress and 
displacement constraints (see eq. (8a) or (8b)). These 
singularities will be referred to as local singularities, since the 
relations are derived by consideration of local segments of the 
structure. This localization is reflected in the banded nature 
of the equilibrium and compatibility matrices, \B] and [C], 
respectively. 

Constraint Dependence for a 20-Bay Truss 

For a more comprehensive example of how easily these 
more subtle forms of singularities can arise, consider the 
20-bay truss shown in figure 2. The truss( 101 ,80) consists 
of 20 bays, each defined by six adjunct elements. The truss 
has 101 elements and has 101 stress and 80 displacement 
degrees of freedom. It has 21 compatibility conditions 
( n — m — 21). Using the evidence found earlier of possible 
global singularities, one should restrict the set of active 
constraints to be no more than the number of displacement 
degrees of freedom (m — 80). The first two cases in table II 
show the rank of the resulting gradient matrix w hen more than 
80 constraints are active. One might consider the restriction 
of limiting the number of active constraints to 80 (out of the 
maximum of 181 prescribed behavior constraints) to be rather 
mild since it might be expected that fewer than 80 constraints 
would be active, anyway. However, this restriction is insuffi- 
cient to prevent all singularities that can arise. 

Consider, first, the stress-displacement relations of the 
structure, which are defined through the modified equilibrium 
matrix [B\. A few typical stress-displacement relations for the 
truss (101,80) are as follows: 

a, = b\ 2 X 2 (9a) 

04 — £ 4 . 5^5 + ^ 4 . 6^6 (9b) 

^88 = ^ 88 , 69^69 + ^ 88 . 70^70 + ^ 88 , 71 ^ 7 ! + ^ 88 . 72^72 ( 9 c) 

<7 9 0 = ^90.69^69 + ^90.73^73 (9d) 

From Figure 2, we can see that the stress in element 1, cq, 
and the vertical displacement at node 2, X 2 , are directly 
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(b) (c) 

(a) Overall view, 

(b) View showing the flexural load path under gravity load. 

(c) View of bay 2 . 

Figure 2— 20-bay truss. Displacement constraint numbering for all nodes t such that 2 < i < 40 is as follows: ( 1 ) for horizontal displacements, the constraint 
index is 2(i - I)- I; (2) for vertical displacements, the constraint index is 2(i - I). (Elements are circled, nodes are not.) 


related as in equation (9a). Similarly, equation (9c) can be 
derived by examining the connectivity of element 88 in the 
18th bay in figure 2. The stress in that element a 88 can be 
determined from the horizontal and vertical components of the 
displacements at node 36, X 69 and X 70 , and at node 37, X 1] 
and X 12 . The other relations can be derived in a similar way. 
Since each relation is derived by examining a small subset of 
the structure, the bandwidth of the f B) matrix will be small, 
and the singularities produced can be considered to be local 
singularities. 

The element areas, which are the design variables for this 
problem, do not appear explicitly in equations (9), since the 
coefficients are independent of those areas. Taking, as 
examples, the sets of constraints {a h X 2 J, {a 4 ,X 5 ,X 6 J, [a 88 ,X 6 9, 
X-jq , X-j \ , Xj 2 1 * and as active constraint sets, would 

produce rank-deficient gradient matrices [Vg] and singular [H] 
matrices (see cases 3 to 6 in table II for numerical veri- 
fication). These clearly demonstrate that subsets containing 
many fewer active constraints than the global singularity 
restriction of at most 80 stress and displacement constraints 
can cause singularities. Next, consider the stress compatibility 


conditions (CC) of the truss(101 ,80), which includes one 
external CC and 20 individual bay CC. The external CC, 
which has 20 terms, has the following form: 

C\' 2°2 **1,7^7 T C {a2 <*\2 »•••* T C | 97(797 = 0 (10a) 

The 20 individual bay compatibility conditions, each of 
which involves six stresses, have the following cyclic form: 

f' 2 , 1^1 0 . 2^2 + C2 3 CJ 3 -T C 2 4O4 + ("2 5(75 + ('2.6^6 = ^ 

(10b. 1) 

Cv6°6 + (*3 7(77 + C 3 8 (J 8 + £39(79 + C3 | 0 ( 7 jq + ("3 j|( 7 n -= 0 

(10b. 2) 

c 2 1 . 96^96 + c 2 ! . 97^97 + * 21 , 98^98 + ? 2 1 , 99^99 + ^ 21 . 100 a 100 

+ ^21.101^101 = 0 (10b. 20) 
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TABLE II.^RANKS OF SEVERAL CONSTRAINT-GRADIENT MATRICES FOR 
THE 20-BAY TRUSS 

|The element areas are pseudorandom numbers between 1.0 and 11.0 in - 1 


Case 

Constraint set 

Rank of 
constraint- 
gradient 
matrix 

Nature of 
singularity 

Number 

Constraints 

1 

181 

] All stresses and displacements) 

80 

Global 

2 

101 

[All stresses) 

80 

Global 

3 

2 

j a | , AT 2 j 

1 

Local 

4 

3 


2 



5 

5 

) a K8.^69’^70'^71'^7:! 

4 



6 

3 


2 



7 

6 

1(7, , 0 6 j 

5 



8 

8 

a 6’ ff 20'^W)' 

7 



9 

6 

'< 

5 



10 

20 

[ .0 1 2 ’ • • ■ '^ 97 ! 

19 
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Again, these relations (eq. (10)) are independent of the 
design variables. The sparsity, and (in the case of the 20 bay 
compatibility conditions) the bandedness, of the [C] matrix 
is produced from the localized character of the relations. The 
external CC (eq. (10a)) ensures that the total deformation 
between the two boundaries (nodes 1 and 41) is zero (see tig. 
2(b)). The 20 bay CC’s ensure that the bars which form the 
bay deform in a consistent manner, such that all bars fit 
together before and after deformations, without inducing any 
residual strains (see, for example, fig. 2(c)). 

Taking the first six stress constraints { a^a 2 o 6 ] will 

produce a constraint-gradient matrix with a rank less than six 
(see case 7, table II). In fact, taking any set of constraints that 
include these six will produce a rank-deficient constraint- 
gradient matrix (e.g., see case 8, table II). Other examples 
of singularities can be produced with sets of constraints that 
include •• or jtf 2 ' a 7 ' <7 i 2 ^--’ a 97 i* These have also 

been verified numerically (see cases 9 and 10, table II). Note, 
here, that the 20 stresses S<t 2 ^ 7 ^i 2 - • correspond to the 
flexural load path for the structure under gravity loads and 
therefore are very likely to all become active, creating 
singularities during structural optimization procedures. 


Discussion 

We have seen how singularities can (sometimes quite easily) 
arise in structural optimization because of rank deficiencies 
in constraint-gradient matrices that are used in mathematical 
programming techniques to determine new search directions. 
In some variants of these algorithms, starting from an identity 
matrix, the matrix [H] or related matrices are constructed by 
using constraint gradients following an update procedure 
(ref. 4). Often these matrices are initially well behaved, but 
become poorly conditioned (or singular) as the solution 
approaches the optimal point, thus requiring re-initialization. 
Noting that, typically, few' constraints are active initially and 


that many may be active near the optimal solution, it is 
reasonable to hypothesize that the singularities observed are 
due to linear dependence among the active constraints. Once 
again, it seems very important to ensure that the active 
constraints chosen constitute a linearly independent set of 
functions. 

Although generating conditions in which singularities can 
arise is very easy using equations (8), producing an algorithm 
to ensure linear functional independence among stress and 
displacement constraints is not so straightforward. One naive 
approach, which would be applicable to any structural analysis 
formulation, might be to find the rank of the matrix formed 
by the gradients of the active constraint set. If this matrix were 
rank deficient, the offending constraint(s) could be deleted, 
or replaced with some “less active, but (what one would hope 
are) independent, constraint(s). 

A far less costly, and more elegant, approach would be to 
examine the stress-displacement relations, the displacement- 
stress relations, and the compatibility conditions of the IFM 
(eqs. (8)). Although the stress-displacement relations alone 
would be sufficient, it would often be advantageous to use one, 
or both, of the other two relationships. For example, when 
only stress constraints occur in the active constraint set, 
examination of the CC alone would be sufficient to determine 
whether or not the set was linearly functionally independent. 
Details on a procedure to ensure linear functional independence 
among active constraints will not be given here, but will be 
provided in a subsequent paper (Guptill, J.D. ; Patnaik, S.N. ; 
and Berke, L.: Identification and Preclusion of Singularities 
during Structural Optimization; to be published). 

Multiple Load Conditions 

With multiple load conditions that are mutually independent, 
it can be shown that the stress gradients for any single element 
are mutually independent. To illustrate this result with a simple 
numerical example, we took two load conditions for the three- 
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TABLE ITI.-RANKS OF SEVERAL CONSTRAINT- 
GRADIENT MATRICES FOR MULTIPLE 
LOAD CONDITIONS 
[Element areas: LOO, 9.30, and 3.67 in. 2 ] 


Case 

Constrain! set 

Rank of 
constraint- 
gradient 
matrix 

Number 

Load 

conditions^ 

1 2 

1 

2 

K! loii 

2 

2 

2 

j (72 1 i X 2 1 

2 

3 

3 

Iffj.ffil r<7.0 

3 


‘H oad conditions arc defined in figure I 


bar truss shown in figure 1. The linear independence of the 
stress gradients associated with the first element under the two 
load conditions is demonstrated in table III. Furthermore, 
although a 2 and X 2 were found to be linearly functionally 
dependent under one load condition, taking a 2 from load 
condition 1 and X 2 from load condition 2 produces a gradient 
matrix of full rank (see table III). Moreover, taking two 
linearly functionally independent constraints from one load 
condition and one from another produces a gradient matrix 
of rank three. An example of this is also shown in table III. 

Frequency Constraints 

Frequency constraints appear to be linearly functionally 
independent of stress and displacement constraints. The fre- 
quency constraint/is inversely proportional to the frequency 
o) where 

2 xliwL fL\s]K 

03 = =«- = -=*= =“=*“ ( 11 ) 

xlmx* Fl\M 7 \F„ 

and X u , [tf], and [M] are the displacement mode shape, the 
stiffness matrix, and the mass matrix, respectively, associated 
with the displacement method (ref. 8); and F^, [5], and [M r \ 
are the force mode shape, the IFM governing matrix, and the 
IFM mass matrix, respectively (ref. 9). It seems unlikely that 
the frequency constraint could ever be written as a linear 
combination of the stress and displacement constraints. 
Numerical examples indicating the linear independence of 
frequency constraint gradients are shown in Table IV. 

Structures Other Than Trusses 

For nontruss structures, the complexity is increased for 
several reasons, including (1) the presence of more than one 
design variable per element, (2) the presence of multi-axial 


TABLE IV. -RANKS OF SEVERAL CONSTRAINT-GRADIENT 
MATRICES IN THE PRESENCE OF FREQUENCY CONSTRAINTS 
[Element areas: 1. 00, 9.30, and 3.67 in. 2 ] 


Case 

Constraint set 

Rank of constraint- 




gradient matrix 


Number 

Constraints 

With 

Without 




frequency 

frequency 




constraints 

constraints 

1 

6 

j(7 | yj ij 

3 

2 

2 

4 

j(7|,02.CTi/! 

3 

2 

3 

3 

1*2 <* 2 / ! 

2 

1 


stress states, and (3) more complicated IFM equations. An 
initial introduction to the complexity of singularity issues for 
nontruss structures has been reported (ref. 9). 


Conclusions 

Singularity conditions can arise in structural optimization 
because of linear functional dependence among active stress 
and displacement constraints. These conditions can be global 
or local in nature. Local singularities can occur more fre- 
quently than global singularities. Linear functional dependence 
can be seen among sets of constraints containing very small 
percentages of the prescribed behavior constraints. 

The presence of linear functional dependence can best be 
determined by examination of the equations derived in the 
integrated force method. If the active constraint set is 
composed of stress constraints only, examination of the IFM 
compatibility conditions is sufficient. When the active con- 
straint set includes both stress and displacement constraints, 
either the stress-displacement relations or a combination of 
displacement-stress relations and the compatibility conditions 
need to be examined. Frequency constraints appear to be 
linearly functionally independent of the stress and displacement 
constraints. Stress and displacement constraints that cross 
multiple (mutually independent) load conditions are linearly 
functionally independent. Although more complicated, it would 
be of benefit to extend this analysis of singularities to nontruss 
structures. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 12, 1991 
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Appendix A 
Symbols 


A, cross-sectional area of the / th element 

[B\ equilibrium matrix 

[B\ modified equilibrium matrix 

h, , (/J)th element of the modified equilibrium matrix 

[C] compatibility matrix 

[C] modified compatibility matrix 

Cjj (ij ) th element of the modified compatibility matrix 

E Young's modulus 

E ) Young's modulus of the /th element 

[£ t ] diagonal matrix of Young’s moduli scaled with 
lengths 

F internal forces 

force mode shape associated with IFM frequency 
analysis 

/ frequency constraint 

[Gj concatenated flexibility matrix 

g } y'th constraint value 

im coefficient matrix in a method of feasible 

directions 

[/] identity matrix 

[/] deformation coefficient matrix 

[7] modified deformation coefficient matrix 

m stiffness matrix associated with the displacement 
method 

f characteristic dimension of the three-bar truss 

i) length of the /th element 

[M\ mass matrix associated with the displacement 
method 

[M t \ IFM mass matrix 

m number of displacement degrees of freedom 

n number of design variables; number of stress 

degrees of freedom 

P load vector 


P* modified load vector 

r number of strain compatibility conditions 

[S] IFM governing matrix 

X displacement vector 

Xj yth nodal displacement 

X u displacement mode shape associated w ith the 

displacement method 

{ctj] set of constants 

(3 deformation vector 

j3 0 initial deformation vector 

5R effective initial deformation vector 

l node number 

k number of active constraints used in the 

calculation of the direction vector 

Xq scaling factor 

a stress vector 

a 0 strength 

<7, stress of the /th element 

X design variable vector 

X, /th design variable 

ip direction vector 

frequency 

0 null vector 

Vf gradient of the merit function 

[V#] constraint-gradient matrix 

Vgj gradient of the yth constraint 

dgj/d\i partial derivative of the /th constraint with 
respect to the /th design variable 

['] T transpose of a matrix 

[ ] _l inverse of a matrix 

[ • ) 7 inverse transpose of a matrix 
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Appendix B 

Key Integrated Force Method Equations 


The key equations of the integrated force method (IFM) as 
related to the singularity issue are presented in this appendix. 
The IFM considers all the internal n forces, F, as the simultan- 
eous unknowns. The m force equilibrium equations, JZ?]/ 7 = P, 
and the r strain compatibility conditions, [C][G]F = 0, are 
concatenated to obtain the governing equations of the method 
(refs. 10 and 1 1) as 


[tf] 

F = 

~ p- 

|C][C]_ 


8R 


or [S]F = P* 


(Bl) 


where [#] is the m x n equilibrium matrix, [C] is the r x n 
compatibility matrix, [G] is the njx n concatenated flexibility 
matrix that links deformations 0 to forces F as (/ 3 = [G]/ 7 ), 
P is the w-component load vector, 5R is the r-component 
effective initial deformation vector defined as 5R = -[C]0o, 
0o is the ^-component initial deformation vector, [S] is the 
n X n governing matrix, and m + r = n. The matrices ffl], 
[C], [G], and [S] are banded and have full row ranks of m, 
r , n, and n , respectively. For simplicity, initial deformations 
are neglected here (0 O = 0). 

The solution of equation (Bl) yields the n forces, F. The 
m displacements, X, are obtained from the forces by 

*=l /IIGIF (B2) 

where [7] is the m X n deformation coefficient matrix defined 
as ([7] = m rows of [S] ~ T ). For static analysis, the matrix [7] 
is independent of element areas. For trusses, the flexibility 
matrix [G] is a diagonal matrix, and its elements are 

£ 

g(i,i)=— — for/ =1,2 n (B3) 

A.E. 


where f, is the length of the ith element, A, is the area of the 
/th element, and E, is the Young’s modulus of the ith element. 

Stress-Displacement Relations 

The stress-displacement relations can be obtained from the 
general relations (0 = [G]F = [5] 7 X) with appropriate spe- 
cialization for trusses as 


a - [B]X (B4) 

where [B] is a sparse, banded matrix with [ B ] = [£ ( j[£] 7 , and 
Ejftj are the nonzero elements of the diagonal matrix [£)], 

Displacement Stress Relations 

The displacements X can be written in terms of the stresses as 

X= [J\o (B5) 


where |7] = [7][£ f ] 

Compatibility Conditions in Terms of Stresses 

The strain compatibility conditions ([CHG]/ 7 ^ 0) can be 
written as stress compatibility conditions: 

[C\o = 0 (B6) 

where [C] is a sparse matrix with [C] = [C][E { ] _l . The [ 5 ], 
[7|, and [C] matrices can all be considered independent of 
element areas for static analysis. 
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